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Kink oscillations of flowing threads in solar prominences 
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ABSTRACT 



Context. Recent observations by Hinode/SOT show that MHD waves and mass flows are simultaneously present in the fine structure 
of solar prominences. 

Aims. We investigate standing kink magnetohydrodynamic (MHD) waves in flowing prominence threads from a theoretical point of 
view. We model a prominence fine structure as a cylindrical magnetic tube embedded in the solar corona with its ends line-tied in 
the photosphere. The magnetic cylinder is composed of a region with dense prominence plasma, which is flowing along the magnetic 
tube, whereas the rest of the flux tube is occupied by coronal plasma. 

Methods. We use the WKB approximation to obtain analytical expressions for the period and the amplitude of the fundamental mode 
as functions of the flow velocity. In addition, we solve the full problem numerically by means of time-dependent simulations. 
Results. We find that both the period and the amplitude of the standing MHD waves vary in time as the prominence thread flows 
along the magnetic structure. The fundamental kink mode is a good description for the time-dependent evolution of the oscillations, 
and the analytical expressions in the WKB approximation are in agreement with the full numerical results. 

Conclusions. The presence of flow modifies the period of the oscillations with respect to the static case. However, for realistic 
flow velocities this effect might fall within the error bars of the observations. The variation of the amplitude due to the flow leads to 
apparent damping or amplification of the oscillations, which could modify the real rate of attenuation caused by an additional damping 
mechanism. 
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1. Introduction 

Recent observational evidence of ubiquito us periodically vary 
ing featur es in the solar corona (e.g., Tomczyk et alj 200 
iJess et alJ[2009t iTomczvk & McIntoslJl2009t IWang ~ 
has raised the debate on whether these observations are caused 
by magnetohyd rodynamic (MHD) waves or b y quasi-periodic 
flows (see, e.g.. iDe Pontieu & Mcinto sh 2010). There seem to 
be strong theoretical arguments supporting the wave interpreta- 
tion (e.g..lErdelvi & Fedunll2007t | Van Doorsselaere et al.l 2008; 
iTerradas et al.l 120101 : IVerth et all l2010h . However, waves and 
flows are not mutually exclusive and, in fact, both phenomena 
have been simultan eously observed in th e fine structure of solar 
prominences (e.g., Okamoto et al. 2007). This offers us the op- 
portunity to study the interaction between waves and flows in the 
solar atmosphere. 

i 

The fine structure of solar prominences is clearly visi- 
ble in the high-resolution Ha and Ca II H-line images from 
the Solar Opti cal Telescope (SOT) aboard th e Hinode sa tellite 
(e.g..lOkamoto et al.1120071: iBerger ■ et dj|2008t IChae et alj |2008t 
iNing et alJl2009t ISchmieder et alJl2010t IChaefeoiOh . When ob- 
served above the limb, vertical structures are commonly seen 
in quiescent pr ominences (e.g., iBerger et alj 120081 IChae et al] 
2008: IChad20 1 Oh . while horizontal threadlike str uctures are usu- 
ally o bserved in active region prominences (e.g., Okamoto et al. 
2007). Although it is ap parently difficult to rec oncile both pic- 
tures, some authors (e.g., Schmie der et alj |2010) have suggested 
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that vertical threads might actually be a pile up of horizontal 
threads which appear as vertical structures when projected on the 
plane of the sky. This idea is consistent with Ha observations of 
filam ents on the solar disk from the Swedish Solar Telescope 
(e.g.. iLin etaill2008l 120091) . in which the filament fine struc- 
ture is seen as thi n and long d ark ribbons. On the other hand, 
other authors (e.g., Chae 2010) have argued that vertical threads 
are real and are an indication of the existence of vertical mag- 
netic fields in quiescent prominences. Thus, it remains unclear 
whether all prominences have the same magnetic sttucture or, 
on the contrary, the magnetic field in quiescent prominences is 
predominantly vertical and active region prominences have hor- 
izontal fields. A recent review on the properties of prominence 
threads can be found in|ni3<|2Q10|). 

There are many evidences of transverse oscillations of 
the fine structures of both active region and quiescent promi- 
nences, which have been interpr eted in terms of kink MHD 
waves (see the recent re views bv !BaUesteill200a lOliverJ 120091 
Arregui & Bal lestetll2010l) . The reported periods are usually in 
a narrow band between 2 and 10 minutes, while the oscilla- 
tions are typically damped after a few periods. In addition, flows 
and mass motions i n prominence s have been also rep o rted (e.g., 
Zirk er etalJ [19981: IWaiigl [19991 iKucera et alJ 120031: ILin et al.l 
120031 : lAhn et al.ll2010l) . The typical flow velocities are less than 
30 km s in quiescent prominences, although larger values up to 
40-50 km s _1 have been observed in active region prominences. 

The work of lOkamoto et al.l (120071) is an example of simul- 
taneous transverse oscillations and mass flows in prominence 
fine structures observed with Hinode/SOT. In the present pa- 
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pe r we focus on the the oretical analysis of the event reported 
bv lOkamoto et al.l {2007). Similar observations of simultaneous 
flows and oscillations have been repor t ed by lOfman & W ang 
(2008) in coronal loops and bv lCao et al.l(l2010[) in filament foot- 
points. Also, the recent work by lAntolin & Verwichtd d201 lb 
on observations of transverse oscillations of loops with coro- 
nal rain is rel e vant f or our present theoretical investigation. 
Oka moto et al.l d2007l) observed an active region prominence 
formed by a myriad of horizontal magnetic flux tubes which 
are partially outlined by threads of cool and dense prominence 
plasma. The magnetic tubes are probably rooted in the solar 
photosphere. Although only the part of the tubes filled with 
prominence material can be seen in the Ca II H-line images, the 
length of the whole magnetic tube must be much longer than 
the length of the prominen ce threads, which is ro ughly between 
3,000 km and 16,000 km. lOkamoto et all d2007l) detected that 
some threads were flowing along the magnetic tubes and simul- 
taneously oscillating in the vertical direction. The mean period 
of the oscillations was 3 min and the apparent flow velocity on 
the plane of sky was around 40 km s _1 . The oscillations were in 
phase along the whole length of the threads, and the wavelength 
was estimated to be at least 250,000 km. 

The event observed by lOkamoto et aL (120071) was stud - 
ied from a theoretical point of view by Terradas et al. (2008), 
who interpreted the observations in terms of standing kink 
MH D modes supported by the magnet i c structure (see, e.g. , 
Edwin & Robertsl ll98l iDfaz et al.ll2002L iGoossens et all 120091 
Soler et all 120101). A n interpretation of the observations by 
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(12001) in terms of kin k modes was also sug- 



gested by lErdelvi & Fedur] (120071) and IVan Doorsselaere et al.l 
(2008). Terrad as et al.l (120081) used the observed wave properties 
provided by |Okamoto et al. (2007) to perform a seismological 
estimation of a lower bound of the prominence Alfven speed. 
The ti me-dependent numerical simulations by iTerradas et al.l 
(2008) suggested that the influence of the flow on the period 
was small. Nevertheless, the precise effect of the flow was not 
assessed in their work because a detailed parametric study was 
not performed. The purpose of this paper is to advance the anal- 
ysis of the event observed by Okam oto et al.l {2007) by com- 
bining both analytical and numerical methods. In the analyti- 
cal part, we use the WKB approximation to assess the effect 
of the flow on the period and the amplitude of the transverse 
oscillations. Expressions of these quantities as functions of the 
relevant parameters of the model are obtained. In the numerical 
part, we go beyond the WKB approximation and solve the full 
time-dependent problem. The implications of our results for the 
magneto-seismology of prominences are also discussed. 

This paper is organized as follows. Section [2] contains the 
model configuration and the basic governing equations. The an- 
alytical investigation of standing kink MHD waves in flowing 
prominence threads using the WKB approximation is included 
in Section [3] while the full numerical solution of the time- 
dependent problem is performed in Section|4] Finally, our results 
are discussed in Section[5] 

2. Model and basic equations 

The background model in which the waves are superimposed is 
schematically shown in Figure Q] It is composed of a straight 
and cylindrical magnetic flux tube of radius R and length L, 
whose ends are fixed at two rigid walls representing line-tying 
at the solar photosphere. The z-axis is chosen so that it coincides 
with the axis of the tube, and the photospheric walls are located 
at z — +L/2. The magnetic tube is partially filled with promi- 



nence plasma of density p p , while the rest of the tube, i.e., the 
evacuated part, is occupied by less dense plasma of density p e . 
The density of the external plasma is the density of the coronal 
medium, p c . The length of the prominence region (thread) is L p . 
The thread flows along the tube as a block with constant speed 
vo- The magnetic field is B = Be,, with B homogeneous. As the 
B — approximation is used in the present work, with B the 
ratio of the gas pressure to the magnetic pressure, the plasma 
temperature is irrelevant for the study of kink MHD waves 
supported by the model. In the absence of flow, standing kink 
M HD waves supported by t he present model were investigated 
by iJoarder et all d!997l) andlDfaz et all (|2001|) in Cartesian ge - 
ometry, and bvlDfaz et al.l (f2002l).lD vmova & Ruderman (2005), 
IDfaz et all (1201 Ol) . and lSoler et al.l d2010l) in cylindrical geome- 
try. 

We adopt the TT approximation, which is valid for R/L <sz 1 
and R/Lp <K 1. To check whether or not this approximation is 
reasonable in the context of prominence threads, we take into 
account that the values of R and Ln reported by the observations 
(e.g., lLinl[2004t lOkamoto et aDl2007l IUnetal] l2008) are in the 
ranges 50 km < R < 300 km and 3,000 km < L p < 28,000 km, 
respectively, and assume L ~ 10 5 km as a typical length for 
the magnetic tube. We obtain R/L p and R/L in the ranges 
3 x 10~ 3 < R/L p < 0.1 and 5 x 10~ 4 < R/L < 3 x 10 -3 , 
meaning that the use of the TT approximation is justified in 
prominence fine structures. In the case vo = 0, the basic equa- 
tion governing linear kink MHD waves of the flux tube in the 
TT approximation was derived bv lDvmova & Rudermanl (f2005) 
in their Equation (21 ). In the absence of flow, the TT approx- 
imatio n wa s used bv iDymova & R udermarj (120051) . [Diaz et al. 
(120101) . and|So!eretaL| (|2010|). The results of these works fully 
agree with the gene ral results be yond the T T approxim ation by 
IJoarder et all d 19971) . IDfaz et alJ (feOQll) . and lDfaz et all d2002l) . 

In the prese nce of flow, an intui t ive g eneralization of 
Eq uation (21) oflbvmova& Rudermarj (120051) was performed 
by Terradas et al. (2008) in their Equation (2). Mathematically, 
iMorton & Erdelvil d2010al) also considered the variation of den- 
sity with time and obtained a sim ilar expression in their 
Equation (18). We refer th e reader to ITerradas et al.l (120081) and 
Morton & Erdelyi (2010a) for a detailed derivation of the basic 
equatio n. In the mathematical derivation of IMorton & Erdelvi 
(2010a) it is assumed that the difference of the flow velocity 
between the internal and external plasma is small, i.e., much 
smaller than the Alfven velocity. So, we restrict our present in- 
vestigation to values of the flow velocity that satisfy Vo / v'a p 1 , 
-jL= is the prominence Alfven speed. Assuming 



where 



'Ap 



B — 50 G and p p = 10~ 10 kg rrT 3 as realistic values of the mag- 
netic field strength and density in active region prominences, we 
obtain VAp ~ 446 k m s -1 . Since the flow v elocities on the plane 
of sky estimated by Okamoto et al.l d2007l) are in the interval be- 
tween 15 km _1 to 46 km s _1 (see their Table 1), the restriction 
vo/vAp <K 1 is satisfied for realistic parameters in prominences. 
Thus, the governing equation we study in the present work is 



d 2 v r (z, t) 
dt 2 



-vl(z,t) 



d 2 v r (z, t) 
dz 2 



= 0, 



(1) 



which has to be solved along with the condition of line-tying at 
the photosphere expressed as v r (+L/2, t) = 0, and a given initial 
condition at t = 0. In Equation ([TJ, v r (z, t) is the radial velocity 
perturbation at the tube boundary and Vk(z, t) is the kink speed, 
which in our model is a function of z and t, namely 



Vk(z, f) 



Vk P if \z-Zo - vot\ < Lp/2, 
vice if \z-Zo~v t\ > L p /2, 



(2) 
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Fig. 1. Sketch of the prominence fine structure model adopted in this work. 



where 



Vkp 



2B 2 



Vke 



2B 2 



\yU(p e +p c )' 



(3) 



with /j the magnetic permittivity, and zo corresponds to the po- 
sition of the center of the prominence thread with respect to 
the center of the magnetic tube at t = 0. Note that zo < if 
the thread is initially located on the left-hand side of the cen- 
ter of the tube, whereas zo > if the thread is initially located 
on the right-hand side. We see that the flow does not explic- 
itly appear in Equation (Q]) since it is enclosed in the d efini- 
tion of Vk(z, t) given in Equation (0. Terrad as"et alJ (|2008j) per- 
formed time-dependent simulations and solved Equation (Q]) nu- 
merically. Here, our aim is to solve Equation ([TJ by using both 
analytical and numerical methods. 



3. Analytical investigation: WKB approximation 

We solve Equation ([TJ analytically by using the Wentzel- 
Kramers-Brillouin (W KB) approximation (see, e.g., 
iBender & Orszagl Il978i for details about the method). The 
WKB approximation has been recently applied t o the investi- 
gation of MHD waves in cooling coronal loops dMorton et alJ 
2010t iMorton & Erdelvil |2010atfbh . In particular, the work by 



Morton & Erdelvil (1201 Oal) is especially relevant for the present 
investigation as they studied kink oscillations of coronal loops 
with variable background. 

The WKB approximation is an approximate method to study 
waves in a changing background whose properties are smooth 
functions of space and/or time. In the present application of 
the WKB approximation we assume that the time scale related 
to the waves, e.g., the period, is much shorter than the time 
scale related to the changes of the background configuration. 
Under these conditions, it is possible to define a time-dependent 
frequency which slowly varies because of the changing back- 
ground. To apply the WKB approximation we define the param- 
eter 5 as 



6=—. 
L 



(4) 



The validity of the WKB approximation is restricted to small 
values of 6 so as P5 <K 1, w here P is the period o f the oscil- 
lations. In the observations by lOkamoto et alJ (|2007), the mean 
flow velocity and period are vo ~ 40 km s _1 and P » 3 min. For 
L ~ 10 s km these values result P5 ~ 0.072, meaning that the 



condition of applicability of the WKB approximation is fulfilled 
in the case of transverse oscillations of flowing threads. 

Using the parameter 5 we define the dimensionless time, tt, 

as 

h = 6t, (5) 
and we express the solution to Equation ([JJ) in the following form 



v r (z,h) = £Mz.*i)exp(~ni(fi)J, 



(6) 



with Q\(z,ti) and Qi(fi) functions to be determined. Next, we 
combine Equations ([JJ) and (0, and separate the different terms 
according to their order with respect to 6. As 6 is small, the dom- 
inant terms are those with the lowest order in 5. We obtain two 
equations for Q\{z,t\) and £l\{t\) taking the terms with 0(d°^ 

and 0(6^, namely 



dz 2 

dt 2 



du 



l(z,h) 



+ 2 



dQ x {z,h) d£i x {t x ) 



du 



du 



0. 



(7) 



(8) 



E quations © and I© are equivalent to Equations (24) and (25) 
of IMorton & Erdelvil d2010al) . respectively. 

Now, we define the time-dependent frequency, co (t\ ), as 



which allows us to rewrite Equation ([7]) as follows 



d 2 Q l (z,t l ) oj 2 (h) 



dz 2 



v 2 k (z,h) 



Qi(z,h) = 0. 



(9) 



(10) 



Equation ([T0l has to be solved taking into account the boundary 
conditions Q\(±Lj2,t\) = due to photospheric line-tying. By 
solving Equation fllOl . the dependence on z of function Qi(z, h) 
can be obtained. In addition, since vi(z, h) is a piecewise con- 
stant function of z (see Equation (O), the analytical solutions to 
Equation ([T0l are trigonometric functions with time-dependent 
arguments. Thus, the general solution to Equation ( fTOb satisfying 
the boundary conditions at z = ±L/2 is 



A 2 (H) cos (^z + 0O-i)) if z- < z < z + , 
A 3 (fi)sin(s£i( Z -§)) ifz> Z+ 
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Fig. 2. Dimensionless frequency, wta p , versus time in units of 
the internal Alfven travel time, ta p = L/va p . Results correspond- 
ing to the fundamental mode and the lowest seven harmonics ob- 
tained by numerically solving Equation (fT4l) for a flowing thread 
with Lp/L = 0.1, v /v Ap = 0.05, and z /L = -0.25. The verti- 
cal dotted line denotes the time when the prominence thread is 
centered within the flux tube. 



of time in Figure [2] for a particular set of parameters. We find 
that the dispersion diagram is symmetric with the time when the 
thread is located at the center of the magnetic tube (denoted by 
a vertical dotted line in Fig. |2]i as point of symmetry. The fun- 
damental mode and the first harmonic are smooth functions of 
time. The other harmonics displayed in Figure [2] show a com- 
plicated set of couplings and avoided crossings. The reason for 
this behavior is that the fundamental mode and the first har- 
monic correspond to global oscillations of the flux tube because 
both the prominence and the evacuated parts of the tube are dis- 
turbed. On the contrary, high harmonics correspond to modes 
more confined within one of these regions. Thus, the collec- 
tio n of modes and th e ir pro perti es are similar to th ose studied 
bv lJoarder & Roberts! d 19921) and lOliver etal] (Q993) in slab ge- 
ometry. 

3.1. Approximation to the fundamental mode frequency 

From hereon we restrict our analysis to the fundamental mode 
of oscillation, whose frequency is the lowest order solution to 
Equation (I14t . To obtain an approximation to the frequency, we 
perform a Taylor expansion of Equation (fT4l) and neglect terms 
with O (lo 4 ^ and higher orders in u. The following expression is 
obtained 



where Ai(fi), A2(fi), and Aj,{t\) are time-dependent coefficients, 
(p{t\) is a time-dependent phase, and z- and z+ denote the loca- 
tions of the interfaces between the prominence thread and the 
evacuated regions, namely 



z±=zq± — + t x L. 

The locations of the interfaces change as the dense thread moves 
along the magentic tube. Q\ must satisfy appropriate boundary 
conditions at z - z±. Since the interfaces corresp ond to con- 
tact discontinuities (see Goedbloe d & Poedtsll2004h . the bound- 
ary conditions are 



[[&]] =0, 



dQx 



dz 



= 0, 



(13) 



where [[X]] stands for the jump of the quantity X at z — z±. 

Applying the conditions of Equation ( fT3l ) on the solutions 
given by Equation ( fTTT ). we arrive at the following equation 



Vkp 



tan 



Vke 



L~ Lr, 



Zt) 



+ t\L 



v kp + v ke cot (sff Lp) tan [g (z + g + t lL )] 
Vkp cot (^Lp) - v ke tan [ig2 ( Zo + Zft + h L)] ' 



(14) 



Equation ( TPfl i is the time-dependent dispersion relation. For 
fixed t\, the solution of Equation ( TPfl i is w(fi). Note that, al- 
though Equation ([Pil l is written in a more compact form, it is 
consistent with dispersion relations previously obtained for nor- 
mal modes in the static case, i.e., vo = 0. Equati on (TT4t with 
vo = is equivalent to Equation (11) of ISoler et al.l (|2010) if the 
substitutions — > — zq and L~ — > + zo are performed 
in their expressi on. Also for vn - , Equation ([Pil l is similar to 
Equation (17) oflJoarder & Roberts! d 19921) and Equation (A5) of 
lOliver et al obtained in Cartesian geometry. 

We have solved Equation (TBI) by standard numerical tech- 
niques. The frequencies of the fundamental mode and of the low- 
est seven harmonics with respect to z are displayed as functions 



oj{h) * 2v kp J — 



(12) yj{L -L P ){L + |L p ) - 4 fa + t x Lf (l + 



(15) 



The effect of the flow is contained in the denominator of 
the right-hand side of Equation ( fTBT ). We see that the effect 
of the flow on the frequency is more complicated than a sim- 
ple Doppler shift. There are two reasons that cause this depen- 
dence. On the one hand, our model is a complicated structure 
in the sense that only the dense prominence material is mov- 
ing. It is well known that a wave propagating in a uniform mag- 
netic tube with a constant siphon flow is affected by a constant 
Doppler shift of the frequency due to the flow. However, the 
effect of the flow is not so simple in more complicated con- 
figurations. Even in the case of a flux tube with a constant 
flow within the tube but no flow in the exterior of the tube the 
wave frequencies suffer correctio ns due to the flow that are not 
simple frequency shifts (s ee, e.g., Nakariak ov & Robertsl [l995: 
iTerra-Homem et al1l2003l) . Our configuration is very different 
from the typical uniform magnetic flux tube with a siphon flow, 
so that the frequency is also modified by the change of position 
of the dense plasma within the magnetic tube. On the other hand, 
we are dealing with standing modes, not propagat ing waves. For 
standing modes in flux tubes iTerradas et al.l (1201 ll) have shown 
that flow produces a spatially dependent phase shift along the 
magnetic tube. In our case, this phase shift is contained in the 
time-dependent phase <p{t\) of Equation (fTTT) . 

For typical prominence and coronal densities, p p » p c 
and p p » p e . Therefore, the term with the ratio of densities 
in the denominator of Equation ( fTBT ) can be neglected. Then, 
Equation (T5[ simplifies to 



(0(h) 



y j(L-L v )(L+±L p )-4(zo + t l L) 2 



(16) 
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In the absence of flow, i.e., t\ = 0, and for zo = 0, Equation (IT6b 
loses its time dependence and becomes, 



CJ : 



2v kp , 



<(L-Lp)(L+ ±L p )l p 



(17) 



Equation (fTTT i is consistent wi th the approximati on of the normal 
mode frequency obtained by iDfaz et aL (2010, Equation (8a)). 
For L p «: L we can approximate L + 



It 

3 ^p 



L, and 



Equati on dTTb reduces to the expression found by ISoler et all 
d2010 . Equation (17)). Although the case L p — > L is very un- 
realistic in prominences because the observed lengths of promi- 
nence threads correspond to L p /L «: 1, it is instructive to take 
into account this limit. For Lp — > L the frequency given in 
Equation ( fT7b tends to infinity. The reason is that the fundamen- 
tal kink mode beha ves as a hybrid mode like those described by 
lOliver et"aT. in a Cartesian slab (se e also the string modes 

investigated bv iJoarder & Roberts! fl992h . Equations (TT3T> — (fT7T> 
ar e approximat i ons o f the hybrid mode frequency. As explained 
bv lOliver et alj d 19931) . hybrid modes owe their existence to the 
presence of both the dense part and the evacuated part of the 
tube. In the limit L p — > L the evacuated part is absent and the 
hybrid kink mode disappears. Thus for L p — > L the fundamen- 
tal kink mode is not the hybrid mode but the first internal mode 
with frequency 



Vkp" 



(18) 



Since L p /L <K 1 in prominences, the fundamental mode in a re- 
alistic thread is the hybrid mode and Equations (fT31> — (TT7b apply. 

We plot in Figure[3]the solution of Equation (TBI) for the fun- 
damental mode as a function of time, along with the analytical 
approximation given by Equation (fTol i. Equation ( TBI has been 
solved by using standard numerical methods to obtain the roots 
of transcendental equations. These computations have been per- 
formed for different values of L p /L and Vo/va p , and for fixed 
Zo/L. There is a very good agreement between the full solution 
(solid lines in Figure[3} and the approximation (symbols). From 
Figure [3 a), we see that the frequency decreases as the length of 
the prominence thread increases. On the other hand, Figure [3jb) 
shows that the variation of the frequency with time is more im- 
portant as the flow velocity gets faster. We find that the mini- 
mum of the frequency takes place when the thread is centered 
within the magnetic tube. Therefore, the minimum of the fre- 
quency depends on both the initial position of the thread and the 
flow velocity since the relation zo + Vot = has to be satisfied. 

Equation dTot corresponds to the instantaneous frequency. 
However, the actual temporal dependence of the oscillation in 
the WKB approximation is given by function Qi(fi). We ob- 
tain Oi(fi) from Equation (O by integrating w(fi) given by 
Equation ( fT6] ). Hence, 



Oi(fi) 



Vkp 



arctan 



2(z +?iD 



- arctan 



A /(L-L p )(L+ \L p )-4(z Q + hL) 

2zo 



^(L-Lp)(L + iL p )-4z2 
where we have used the condition Q.\ — at t\ — 0. 



(19) 



3.2. Dependence of the amplitude on time 

Here, we estimate the variation of the amplitude of the os- 
cillations with time. To do so, we use Equation ([8]). By tak- 
ing into account the definition of the time-dependent frequency 
(Equation (0), we rewrite Equation ([8]l as 



2w(fi) dti 



Qi(z,h) = 0. 



(20) 



Next, we consider the approximate w(?i) for the fundamental 
mode obtained in Equation (fT~5T > to express this last Equation as 



ggifoft) 
du 



2L(z +hL) 



(L-Lp)(L+iLp)-4(z () + fiL) 



:Gifcfi) = 0. 



(21) 



Note that to solve Equation d2Tb we do not have to care about 
the dependence of Q\ on z. For a given z, Equation d2Tb can be 
integrated to obtain the temporal dependence of Q\ at a fixed 
position, namely 



Gifcfi) = Qo(z) 



(L-L p )(L+ iL p )-4(zo + f I L) 2 



(L-Lp)(z,+ iL p )-4z; 



1/4 



(22) 



where Qo(z) is the amplitude at fi = 0. Thus, Equation ( fTTT ) 
gives the spatial dependence of Q\ for a fixed time t\, while 
Equation (l22l provides the temporal dependence at a fixed po- 
sition z. By comparing Equations ( fTol l and d22b , we see that 
Q\(z,t) oc oj(t\Y l l 2 . Then, the amplitude of the oscillation de- 
creases when the thread flows from the center of the tube to the 
footpoint, and increases otherwise. This means that the ampli- 
tude is maximal when the thread is located at the center of the 
magnetic tube. 

We must bear in mind that Equation d22l was derived from 
the terms with <9 in the governing equation, whereas the 

leading terms when 5 is a small parameter are those with (9(<S°). 
Therefore, we have to be cautious about the actual accuracy of 
Equation d22l . although we expect the behavior of the amplitude 
with time to be, at least, qualitatively described by Equation (l22l . 

3.3. Application to magneto-seismology 

The analytical expression of the fundamental mode frequency 
found before can be used to perform magneto-seismology of 
prominence fine structures by using observed periods of oscil- 
lations in flowing threads. We use Equation (fTol l to compute the 
period of the oscillation as a function of time, namely 



Pi(f) 



2n 
u>(t) 



* ^ ^ " [ L + 3^) " 4 (Z ° + V ° ?)2 ' (23) 

which we have explicitly written in terms of the dimensional 
time, f, and the flow velocity, vo. By using Equation ( l23l along 
with observational values of the period, it is possible to give an 
estimation of L, i.e., the total length of the flux tube, which is 
a parameter difficult to measure from the observations. Let us 
assume that we have performed an observation of a transversely 
oscillating and flowing thread with a good cadence and we have 
determined the evolution of the period with time. For conve- 
nience, we set t = when the maximum of the period takes 
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Fig. 3. Dimensionless frequency, ojt\ p , versus time in units of the internal Alfven travel time, ta p = L/va p , for (a) L p /L = 0.05, 
0.1, 0.2 with vo/v Ap = 0.05, and (b) v /v Ap = 0.03, 0.05, 0.1 with L p /L = 0.1. In all cases, zo/L = -0.25. The solid lines are the 
results obtained by numerically solving Equation (1141) . whereas the symbols correspond to the analytical approximation given by 
Equation (1161 1, The vertical dotted line in panel (a) denotes the minimum of the curves, i.e., the time when the prominence thread is 
centered within the flux tube. 



place, namely Pi(0), so we can also fix zo — in Equation ( 1231 . 
Then, we denote as P\(t) the instantaneous period measured at 
t — t. We use Equation d23l and compute the ratio P\(f)l P\(Q) 
to find an estimation of the length of the magnetic flux tube as 



Lp + 



2 2 



1 



(Pl(T)Y 

\PI(P)) 



(24) 



where again we assumed that the flow velocity is slow. We see 
that the right-hand side of Equation d24l > depends on quantities 
that can be directly measured from the observations. As an ex- 
ample, let us assume the following values: Pi(t)/Pi(0) = 0.9, 
vo = 40 km s _1 , L p = 10,000 km, and t — 5 min. From 
Equation (l24l i we obtain L m 1.6 x 10 5 km. However, the ac- 
curacy of Equation d24l i is limited by the uncertainties and error 
bars of the observations. In particular, a very accurate determina- 
tion of the ratio fi(r)/Pi(0) is needed. For instance, a 10% un- 
certainty of Pi(r)/fi(0) produces a 85% uncertainty of L when 
propagation of errors is used in Equation ( 124b and the remaining 
parameters are kept constant. This makes a reliable application 
of Equation d24l) difficult in practice. 

Another relevant parameter that can give us a seismological 
determination of L is the ratio P\j2P 2 , where Pi is the period 
of the first harmonic. The deviation of this ratio from unity is an 
indication of the longitudinal inhomogeneity length scale of the 
magnetic tube. Its application was use d for the first t i me in the 
context of coronal loop oscillations by lAndries et alJ d2005albh 
and has b een explored in su bsequent works (see the recent re- 
view bv lAndries etaLll2009l and reference s therein). In promi- 
nence thread oscillations, Diaz et aD (12010) explored the impor- 
tance of the ratio P\j2P 2 to estimate L p /L in s tatic t hreads. 
While in coronal loops Pi/2P 2 < l. lDfazet all d2010t) found 
that in prominence threads P1/2P2 > 1. The reason for this re- 
sult is that in prominence threads mass density is arranged just 
in the opposite way to that in coronal loops. In loops the density 
is larger in the footpoints than in the apex due to gravitational 
stratification, while in prominence threads the density is much 
larger in the center of the magnetic tube because of the presence 
of the prominence material. 




Fig. 4. Ratio P\/2P 2 of a flowing thread with L p /L - 0.1 and 
zo/L = -0.25 for v /v Ap = 0.03, 0.05, 0.1. The dotted line is 
the result from Equation (f25b in the case without flow and for 
a prominence thread located at the center of the magnetic struc- 
ture. 

For a large density contrast between the prominence and the 
coronal plasmas, and assuming that the thread is located at the 
center of the magne tic cylinder, the relation between P\/2P 2 and 
Lp/L obtained bv lDfaz et al.l d2010h in their Equation (1 1) is 



Pj_ 

2P 2 



(25) 



Let us see how the ratio P] [2Pi K a ffected by the flow and 
so how the results of iDfaz et al.l d2010t) are modified. However, 
it is difficult to obtain an analytical expression for P 2 when 
flow is present. Instead, we compute both Pi and P 2 by solving 
Equation ( fT4l > with numerical methods. Figure [4] shows P\j2P 2 
as a function of time for L p /L = 0.1 and zo/L - -0.25, and for 
different values of vo/va p - The n umerical results a re compared 
with the analytical expression of IDfaz et al.l (120101) . First of all, 
we see that the period ratio is strongly influenced by the veloc- 



R. Soler and M. Goossens: Kink oscillations of flowing prominence threads 



7 



ity at which the prominence thread flows along the flux tube. 
Px/IP-i is maximal when the thread in centered within the tube 
(Pi /2P2 ~ 2.5 for the particular set of parameter in Fig. [4} and 
P1/2P2 — * 1 when the t hread reaches the footpoint. The an- 
alytical approximation of iDfaz eTa □ d2oToh for the static case 
gives a larger value of P1/2P2 in comparison to the case with 
flow. This means that flow reduces the period ratio. Therefore, 
Equations ( l24b and ( l25l ) could be used together to obtain a more 
accurate determination of the magnetic tube length, as the value 
of L inferred from Equation d25l l should be considered as a upper 
bound for this parameter. 



4. Numerical results: time-dependent simulations 

Here, we compare the analytical results of the WKB approxi- 
mation with the full numerical soluti on of the ti me-dependent 
problem. We use the PDE2D code (ISewellll2005h for that pur- 
pose. The set-up of the numerical code is similar to that of 
iTerradas et al.l d2008l) . Equation ([]]) is integrated assuming the 
boundary conditions v r (+L/2, t) = 0. An initial condition for 
v r at t = is also provided. In the code, the distances are ex- 
pressed in units of L and the velocities in units of the prominence 
Alfven speed, va p = ■ We take L — 10 5 km. Assuming 

B = 50 G and p p = 10~ 10 kg rrT 3 as realistic values in active 
region prominences, we obtain va p ~ 44 6 km s -1 . The flow ve- 
locities on the plane of sky estimated by lOkamoto et al.l d2007l) 
are in the interval between 15 km _1 to 46 km s _1 . Hence in our 
simulations we consider values for the ratio vo/Va p in the range 
0.03 < Vo/vAp ^ 0.1. In the code, time is expressed in units of 
the Alfven travel time, i.e., ta p = L/v^p ~ 3.74 min. In addition, 
in all the following computations we have used p p /p c = 200 and 

Pe/Pc = 1. 



4.1. Excitation of the fundamental mode 

First, we use the eigenfunction of the fundamental kink mode as 
the initial condition for v r at t = 0. The eigenfunction is obtained 
by solving the dispersion relation of the normal mode problem 
and computing the spati al distribution of the corresponding per- 
turba t ion (see details inlDvmova & Rud erman 2 0051 ISoler et al.l 
l2010HDiaz etai1l2010h . Hence, we make sure that, after the ini- 
tial excitation, the magnetic tube mainly oscillates in its funda- 
mental mode. 

As a check of the numerical code, we consider the static case 
and put the thread at the center of the magnetic tube, i.e., vo = 
and zo = 0. In this test simulation, we take L p /L = 0.1. By 
looking at the time-dependent evolution of v r , we check that the 
magnetic tube oscillates as a whole. A plot of v r at z = as a 
function of time (not displayed here for the sake of simplicity) 
shows that the amplitude of the oscillation is constant during the 
whole duration of the simulation, meaning that numerical dissi- 
pation is negligible in the simulati on. Later, we perform a power 
spectrum of v r at z — (see, e.g., Carbonell & Ball esterlll99ll) 
and find a large peak centered around the fundamental normal 
mode frequency. The maximum of the power spectrum agrees 
very well with the approximate frequency of the normal mode 
given by Equation ([TBI . For the set of parameters used in this nu- 
merical test, the period in dimensional units is P « 2.6 min. We 
have performed similar simulations but using different values of 
L p /L and zo- Equivalent results to those commented before have 
been obtained in all cases. This indicates, on the one hand, that 
the normal mode interpretation is a very good representation for 



the time-dependent evolution of the oscillation in the static case 
and, on the other hand, that the numerical code works properly. 

Hereafter, we incorporate the effect of the flow. First, we fix 
Lp/L = 0.1 and zo = 0, and consider vq/va p = 0.05. In this situ- 
ation, the thread is initially located at the center of the magnetic 
tube. We make sure that the simulation stops before the threads 
reaches the photospheric wall. We plot in Figure |3a) the radial 
velocity perturbation at z — as a function of time. Contrary 
to the static case, now the amplitude of the oscillation decreases 
in time as the thread moves from the center towards the end of 
the magnetic tube. This behavior is qualitatively described by 
Equation (l22l) obtained in the WKB approximation (see dashed 
line in Fig. [2 a)), although Equation d22l underestimates the ac- 
tual decrease of the amplitude. Inspired by Equation ( l22l . we 
propose the following fit for the amplitude, 



Qi(z,h) = Qo(z) 



(L-Lp)(L+\Lp)-4(z + hL) 2 



{L-Lp){L+\Lp)-Az t 



(26) 



with n an empirical exponent. When Equation ( |26| | is applied 
to the results of Figure |2a), we obtain that the exponent n - 
1 provides a good fit for the amplitude (see the dotted line in 
Figure |3 a)). 

On the other hand, we perform in Fig ure |5fb) a wavelet 
power spectrum dTorrence & Comp"olll998l) of the signal dis- 
played in Figure|5ja). We find that the period decreases in time as 
the thread moves towards the footpoint of the magnetic structure. 
The WKB approximation for the period given by Equation d23l 
is in excellent agreement with the position of the maximum of 
the wavelet spectrum (see dashed line in Figure |5jb)). As dis- 
cussed in Section 13.21 the WKB approximation for the period 
is much more accurate than the WKB approximation for the 
amplitude. In addition, we see that the rate at which the period 
changes with respect to the value at t — is not constant. Using 
Equation d23l and considering the parameters of this particular 
simulation, the period of the oscillation when the thread is at 
z-L/A has decreased of about 9% with the respect to its initial 
value at z — 0, whereas the period decreases of about 45% when 
the thread finally reaches the footpoint of the magnetic tube at 
z = L. 

We repeat the simulation for Lp/L - 0.2 and a faster flow, 
vo/vAp = 0.1, and take zo/L = -0.25 to consider the thread 
initially displaced from the center of the flux tube. The radial 
velocity perturbation at z — is plotted in Figure 0c), and the 
corresponding wavelet power spectrum is shown in Figure |3d). 
These results are equivalent to those of Figures |3a)-(b), i.e., 
both the amplitude and the period of the oscillation depend on 
the position of the prominence thread within the flux tube, taking 
both of them their maximum value when the thread is centered. 
As before, Equation d23l is a very good approximation to the 
period. 



4.2. Arbitrary excitation 

In the previous Section, we have used an initial condition for 
v r that corresponds to the fundamental mode eigenfunction, so 
only this mode is excited. However, it is expected that the en- 
ergy from an arbitrary disturbance of the flux tube is deposited 
in many normal mo des (see a discussion on this issue in, e.g., 
ITerradas et al.i 2007). To represent an arbitrary disturbance of the 
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Fig.5. (a) v r atz = O(solidline in arbitrary units) as a function of time in a prominence fine structure with L p /L = 0.1,vo/va p = 0.05, 
and zo = 0. The initial excitation is the fundamental kink mode eigenfunction for vrj = 0. The dashed line is the amplitude in the 
WKB approximation (Equation (l22l). while the dotted line corresponds to the fit proposed in Equation d26*b with n = 1 . (b) Wavelet 
power spectrum for the dimensionless period, P/ta p , corresponding to the signal displayed in panel (a). The white dashed line is 
the period in the WKB approximation (Equation (l23ll). whereas the horizontal dotted line is the period at t = 0. The red solid line 
denotes 99% of confidence level, (c) Same as panel (a) but for L p /L = 0.2, vo/vap = 0.1, and zo/L = -0.25. (d) Same as panel (b) 
but for the signal of panel (c), and with the horizontal dotted line denoting the maximum of the period. 



flux tube, we consider a Gaussian function as the initial condi- 
tion of v r at t — 0, namely 



v r (z, t = 0) = exp 



(z-0 Z 



(27) 



where f and cr are arbitrary parameters. Whereas f correspond 
to the position of the maximum of the excitation, cr determines 
its width. 

In the following simulations, we consider the same model 
parameters as in Figures |3c)-(d), i.e., L p /L = 0.2, vo/^Ap = 
0.1, and zo/L = -0.25, but use the initial condition given by 
Equation (l27l >. To begin with, we take cr/L — 0.2 and consider 
different values of (. First we use £IL = -0.25, so the excitation 
is mainly confined to the dense prominence region of the flux 
tube. The result of this simulation is displayed in Figure l^a), 
which shows the evolution in time of v r at z = 0, whereas 
Figure |6jb) shows the corresponding wavelet power spectrum. 
It is interesting to compare Figures |5jd) and|6jb) to see that, in 
the present case, the oscillation dynamics is still governed by 
the fundamental normal mode. We see in Figure |6j a) that there 
is some contribution of higher harmonics to the behavior of v r in 
time, although their contribution to the overall oscillation is of 
very minor importance. In addition, the evolution of the ampli- 
tude of v r remains qualitatively described by Equation (l26l > with 
n = l. 

Next, we perform another simulation by taking the same pa- 
rameters as before but assuming g/L = 0. In this case, the max- 
imum of the initial excitation is located in the evacuated part of 



the magnetic tube. Again, we plot in Figure[6]c) v r at z — ver- 
sus time, and in Figure |6jd) the wavelet power spectrum. The 
behavior of v r in time is substantially different in the present 
situation compared to the case of Figures |6ja)-(b). First of all, 
we see that v r is not governed by the fundamental mode ex- 
clusively. The wavelet power spectrum indicates that the energy 
from the initial disturbance is mainly deposited to the fundamen- 
tal mode, but also the first harmonic is excited. The dependence 
in time of the fundamental mode period is again well described 
by Equation d23l >. However, the contribution of the first harmonic 
to the overall oscillation seems to have disappeared when the 
thread is located at the center of the magnetic tube. The reason 
for this result is that the first harmonic eigenfunction has a node 
at z — when the thread is centered, and so the first harmonic 
does not contribute to the signal displayed in Figure |6jc) in such 
a case. On the other hand, it is now difficult to determine the 
effect of the flow on the amplitude of the oscillation. 

Finally, we have performed several simulations for different 
values of cr. If the maximum of the excitation is located within 
the dense part of the flow tube, the results are rather insensitive to 
cr unless values much smaller than L v are used. In all the cases, 
the fundamental mode is predominantly excited. However, the 
results are more affected by the value of cr if the maximum of 
the excitation is located in the evacuated part of the tube. In such 
a case, the larger cr, the more energy is deposited in the funda- 
mental mode. On the contrary, as cr gets smaller the energy of the 
initial excitation is more distributed among higher harmonics. 
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Fig. 6. (a) Same as Figure [2c) but for a initial excitation given by Equation < f27b with g/L = -0.25 and cr/L = 0.2. (b) Wavelet 
power spectrum corresponding to the signal displayed in panel (a), (c) Same as panel (a) but for £/L = and cr/L = 0.2. (d) Same 
as panel (b) but for the signal displayed in panel (c), with the dash-dotted line denoting the period of the first harmonic. 



The results of this Section point out that the time-dependent 
behavior of standing kink MHD waves of flowing prominence 
threads is strongly influenced by the form of the initial distur- 
bance. If the initial disturbance mainly perturbs the dense promi- 
nence part of the flux tube, the oscillations are governed by the 
fundamental kink mode. In such a case, the dependence of both 
the period and the amplitude with the flow velocity are approx- 
imately given by Equations d23l l and (|26| >. respectively. On the 
contrary, the behavior is more complicated if the initial pertur- 
bation takes place in the evacuated part of the fine structure as 
other harmonics are excited in addition to the fundamental mode. 
The contribution of the different harmonics depends on both the 
position and the width of the initial excitation, while the ampli- 
tude of the oscillation does not have a simple dependence on the 
flow velocity. 

5. Discussion and conclusions 

In this paper, we have investigated standing kink MHD waves in 
the fine structure of solar prominences, modeled as coronal mag- 
netic flux tubes partially filled with flowing threads of promi- 
nence material. The present study exten ds and complements the 
previous work by Terradas et al. (2008), who restricted them- 
selves to the numerical investigation of this phenomenon and 
did not perform an in-depth parametric study. Here, we have 
combined analytical methods based on the WKB approximation 
with time-dependent numerical simulations to assess the precise 
effect of the flow on both the period and the amplitude of the 
fundamental kink mode. 

As for the effect of the flow on the period, we can distinguish 
two different situations. On the one hand, we find that the flow 
has a small effect on the period when the thread is located near 



the center of the supporting magnetic flux tube. In this case, the 
variation of the period with respect to the static case may fall 
within the error bars of the observations, and so the effect may 
be und etectable. There our re sults confirm the qualitative discus- 
sion of Terrad as"et al.l (|2008) about the effect of the flow on the 
period. On the other hand, the variation of the period is much 
more important when the thread approaches the footpoint of the 
magnetic structure. Then, the decrease of the period can be larger 
than 50% with respect to the static case. The case in which the 
thread is ne ar one of the footpoint s of the magnetic tube was not 
analyzed by Terrada s~et al.l d2008l) . 

We have also found that the flow affects the amplitude 
of the fundamental mode. This result was not discussed by 
iTerradas et al.l d2008l) . During the motion of the prominence 
thread along the magnetic structure, we find that the ampli- 
tude grows as the thread gets closer to the center of the tube 
and decreases otherwise. This produces an apparent amplifica- 
tion or damping of the oscillations, respectively. Observations 
often indicate that th read transverse oscillations are strongly 
damped (see, e.g.. lLinll200l.l2010tlNin^ et al.l2009l) . While sev- 
eral mechanisms have been proposed and investigated t o ex- 
plain the quick attenuation (se e the recent reviews by lOliverl 
2009: iArregui & Ballest erl2010h . the proces s of resonant absorp- 
tion seems the mo st likely explanation (e.g jArregui et all 2008; 
ISoler et alJl2009allR l2010h . Our present results indicate that the 
actual damping rate of the oscillations might be affected by the 
change of the amplitude due to the flow. This fact should be 
taken into account when the damping rate is used as a seismo- 
logical tool to infer physical parameters of prominence threads, 
because the presence of flow may introduce some uncertaint ies 
on these estimations (see details in lArregui & Ballesterll20ldl) . 
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In addition, our numerical simulations have allowed us to de- 
termine how different perturbations excite the oscillations of the 
magnetic structure. Based on the cases studied in this paper, we 
have obtained that the fundamental mode is mostly excited when 
the perturbation initially disturbs the dense, prominence part of 
the tube. From the wavelet power spectrum of the radial velocity 
perturbation, we conclude that the contribution of higher har- 
monics is negligible, thus the overall oscillation is governed by 
the fundamental mode. On the contrary, a perturbation located at 
the evacuated part of the tube excites the fundamental mode and 
higher harmonics, producing a more complex behavior of the os- 
cillations. In this last case, the effect of the flow on the amplitude 
is more complicated and no simple dependence can be extracted 
from the simulations. 

This paper has explored the properties of MHD waves in 
a coronal magnetic structure with a changing configuration . 
Previous similar work s in this line are , e.g., [T erradas et al. (2008) 
in prom inences, and iMorton et a l. (2010); Morton & Erdelvi 
(201033) in coronal loops. During the revision of this paper 
it also came to our knowledge the recent work by Ruderman 
(1201 11) . In view of the highly dynamic nature of the coronal 
medium in general, and the prominence plasma in particular, this 
kind of modeling represents a better description of the actual 
oscillatory phenomena in the corona and in prominences. The 
present investigation could be extended in the future by incorpo- 
rating the effect of the density inhomogeneity in the transverse 
direction and so investigating the resonant damping of the kink 
mode. 
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